Study of 7-charge correlation in heavy ion collisions, various approaches. 
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Event-by-event 7 — ch correlation is used in studying systems going through QCD chiral phase 
transition. In this paper various methods for measuring 7 — ch correlation in heavy ion collisions 
have been discussed. Dynamical fluctuation due to formation of domains of DCC that can affect 
7 — ch correlation has been discussed. We study known detector and statistical effects involved in 
these measurements and suggest suitable robust observables Avd yn and r m ,i sensitive to small 7 — ch 
correlation signal. These variables are constructed based on moments of multiplicity distributions 
of photon and charged particles. Estimations of expected measurable signals of 7 — ch correlation 
from various available models such as for ideal Boltzmann gas of pions, monte-carlo models based 
on transport and mini-jets have been discussed. Collision centrality dependence of the observables 
have been estimated from Central Limit Theorem and found to be consistent with the model predic- 
tions. We find that observables show high sensitivity to fraction of DCC events and have nonlinear 
dependence on fraction of pions carrying DCC signals. Variation of r m ,l with orders of its higher 
moments m is a observable to extract the nature and strength of 7 — ch correlation. 



I. INTRODUCTION 

Based on decades of experimental searches and theo- 
retical studies it is widely believed that high energy heavy 
ion collisions produce realistic scenario for studying the 
phase transition from hadronic matter to Quark gluon 
plasma. It is believed to be associated with two differ- 
ent transitions, de-confinement and restoration of QCD 
chiral symmetry. Fluctuation of conserved quantities has 
been proposed [l[ to be an important experimental sig- 
nature for such phase transition. Hadronic system that is 
mostly dominated by pions is expected to show a global 
isospin conservation. In such scenario the event-by-event 
isospin number fluctuation is an interesting observable. 
The QCD chiral phase transition is associated with melt- 
ing of 4- vector condensates. An interesting phenom- 
ena like formation of mctastablc domains of "Disoriented 
Chiral Condensate" (DCC) is predicted to occur due to 
the orientation of this condensate relative to the direc- 
tion of its scalar component. Such a phenomena is 
sible for a scenario of rapid cooling like quenching 
for system going from chiral symmetry restored phase to 
broken phase. Formation of DCC domains causes anoma- 
lous production of charged or neutral pions depending on 
the orientation of vacuum towards its pseudo scalar com- 
ponent. Such phenomena might survive final state inter- 
actions and appear in the form of multiplicity fluctuation 
of pions of relative isospinsj3|- As the detected charged 
and neutral particles are mostly from the charged pions 
and the decay of neutral pions respectively this would 
appear in the form of 7 — ch anti-correlation. The ex- 
perimental searches of DCC so far includes searches in 
pp collisions Hi, cosmic ray events and in heavy ion 
collisions |8L fl2 j . Our discussion would be relevant to the 
search in heavy-ion collisions. The prediction for a hot 
medium described by the linear sigma model (l3j showed 
that in the case of central collision of Pb-Pb at SPS ener- 
gies, the likelihood of the DCC events is less than 10 -3 . 
Experimental searches at SPS WA98 experiment 



at y/s=17.3 GeV estimated an upper limit of 3 x ICC 3 . 
It has been argued 14[ that in case of rapid cooling like 
quenching scenario, higher collision energies correspond- 
ing to lower chemical potential (e.g. hrhic < Hsps) 
provides faster cooling rate (\dT/dt\). This suggests that 
RHIC and LHC collisions provide more favorable condi- 
tion for DCC production than SPS collisions. 

From experimental point of view such a study is as- 
sociated with simultaneous measurement of photons and 
charged particles in common phase space with very high 
sensitivity at low momentum. This is because the decay 
of domains of DCC are final stage phenomena of the evo- 
lution of heavy ion collision and the pions carrying signals 
are expected to be of low momentum. A combination of 
pre-shower Photon Multiplicity detectorfPMD)(l5j and 
forward time projection chamber(FTPC) [l6| at STAR 
experiment at RHIC and Photon Multiplicity detec- 
tor(PMD) and Forward Multiplicity Detector(FMD) at 
ALICE experiment [I?) at LHC have the required criteria 
to satisfy such goal. 

In this paper we would like to highlight few issues 
associated to 7 — ch correlation analysis and propose a 
method. We use generating function approach to calcu- 
late different variables and include various detection ef- 
fects like efficiencies, effect of mis-identification etc. Ob- 
servables of 7 — ch correlation are constructed to be suit- 
able for heavy ion collisions that can disentangle dynam- 
ical fluctuation. Assuming formation of DCC domains 
to be one of the probable sources of dynamical signal of 
isospin fluctuation we discuss the sensitivity of the ob- 
servables to the fraction of DCC events and the fraction 
of DCC candidates in an event. Relevant to the heavy 
ion collisions we discuss the centrality dependance of the 
variables. We estimate 7 — ch correlation from various 
models and implement a DCC-modcl based on HIJING 
event generator. 

In the section [IT] we outline the method of construction 
of the observables and their values for DCC events of 
varying fraction. Section IIII1 IIVI and [V] describe the de- 
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tcctor effect like mis-identification, the role of resonances 
and centrality dependance respectively on the proposed 
variables. In sectionlVJwe have calculated the sensitivity 
of the variables on DCC event fraction and pion fractions 
in DCC events. For studying the experimental sensitiv- 
ity of DCC, we have studied various non-DCC models 
in section IVHI and implemented DCC in a Monte-Carlo 
based events in section IVIIII We summarize in section 

ED 



II. METHOD 



Fluctuation of particle ratios has been addressed previ- 
ously in case of conserved quantities like net strangeness 
in terms of kaon-to-pion ratio and net baryons in terms 
of proton to pion ratios. Relevant to our case is the 
study of photon to charge particle multiplicity ratio. Ob- 
servables used in such cases are designed in such a way 
so as to eliminate the statistical fluctuations and at the 
same time be robust against detector inefficiency. A sim- 
ple way of implementing detector efficiencies in terms 
of a binomial probability distribution function say of 
the form P(n,N,e) = N C n e n (l — e) N ~ n would reveal 
the fact that the second moment of observed multiplic- 
ity n is not proportional to second moment of produced 
multiplicity N. The efficiency term e does not factor- 
ize for quantities like variance, skewness and kurtosis . 
However the quantities like observed second and higher 
order factorial moments comes out to be proportional 
to the measured corresponding factorial moments like 
(n(n — 1)) = e 2 (N(N — 1)). Ratios of various factorial 
moments with powers of mean multiplicity would sim- 
ply cancel the explicit efficiency dependence. In case of 
correlation of multiplicities, there could be more compli- 
cated detector effects like mis-identification of one species 
in the form of another, decay and resonance production. 
This could lead to spurious correlation affecting the fi- 
nal results. Also in case of heavy ion collisions there 
are centrality and system size dependence. If heavy-ion 
collisions are assumed to be linear superpositions of mul- 
tiple hadronic collisions, then variables are supposed to 
show number of source scaling (23|. Based on similar 
context and considering various other aspects of par- 
ticle ratio-fluctuation, two observables were introduced 
earlier as measure of 7 — ch-ratio fluctuations. ^dyn^ 
was introduced in Ref [H] and used by STAR Collab- 



oration [19l I20B and 



7— ch 



was introduced by Minimax 



1 variable D = 4 (AQ 2 ) /iVgj, where Q 2 is the variance of the net 
charge(Af + — AT_) in RefQ gives different values for QGP and 
pion gas but depends on efficiency. 



collaboration [2l|. The variable Vdyn is defined as 
7 -ch _ (N ch (N ch - 1)) (A 7 (A 7 - 1)) _ (A ch A 7 ) 



ch 



(iV 7 



w (^ch) 



. (1) 

which for Poissonian case should give zero. The variable 
r m i is defined as 



7 -ch _ (N ch (N ch - 1).. (A ch - m + 1) A 7 ) (A ch ) 



' m,l 



(N ch (N ch - l)..(N ch - m)) (N~ 



(2) 



It is designed such that for all the moments it gives a 
value equal to 1 for Poisson case and higher order mo- 
ments show larger sensitivity to signals. In this section 
we would like to discuss the applicability, robustness and 
sensitivity of these two variable for 7 — ch correlation. 
Since we are interested in fluctuation of ratio of multi- 
plicities let us consider / = N n o/(N n o + N n ±) to be the 
neutral pion fraction. The idea is that by using proper 
combinations of moments we can eliminate the efficiency 
dependence and express our observable in terms of the 
fluctuation of the fraction /. The most efficient way 
of studying the moments including the dynamical and 
detector effect is to follow the generating function ap- 
proach [2ll | where we define, 



G(z) = z " P ( N ) 



(3) 



N=0 



where P{N) denotes the distribution of parent multi- 
plicity where , N = Nq + N c h denotes sum of all neu- 
tral and charged pions. Different moments are calcu- 
lated by taking derivatives of G(z) w.r.to z evaluated at 
z = 1. Considering the fact that the neutral pions are dis- 
tributed according to the probability V(f) the generating 
function has to be modified accordingly 



/v 



1 

G(z ch ,z )= [ dfV(f)J2P(N)[fz + 

N 

(4) 

The distribution V(f) is the event-by-event measured 
distribution of neutral pion fraction. Isospin symmetry 
for a pion gas corresponds to a generic case of pion pro- 
ductions for which V(f) = S(f— 1/3). In case of DCC like 
events0, Q we have V{f) = 1/2 v 7 /- For propagation of 
generating function to include the decay of neutral pions 
to observed photons we apply the "cluster decay theo- 
rem" [22[ . We can express the overall generating function 
as 



Gobs (z c h, z-y) = G (g ch (z c h) , go (z 7 )) 



(5) 



where go(z-y) = z 7 considering the fact that every neutral 
cluster decays into two photons and the charge particles 
do not decay, g c h(z c h) = ^ch- To make the scenario more 
realistic and taking the advantage of same theorem, one 
can include detection efficiencies into the final form of 
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generating function. We consider the observing and non- 
observing as different decay modes with probability equal 
to the detection efficiency So for charged and neutral 
clusters we redefine 



g c h(zch) 



(1 - £ c h) + Ech^ch 

((1 - e 7 ) + e-yZy)- 



(6) 



Here e c h is the efficiency of charge particle detection and 
e 7 is the efficiency of detecting a photon coming from de- 
cay of a neutral pion. We can calculate various factorial 
moments of multiplicity with detector efficiency folded in 
terms of derivatives of final generating function. We can 
define a generalized factorial moment as 



f, 



d m ' n G hs( z ch, z -y) 



dz™ dz™ 



N ch \ NJ. 



-AC^ch- m)! (iV 7 -n)! 

_h — **y — J- 

(7) 

It is convenient to express our observables given in eq[T] 
and eqH]in terms f TOj „ as 



"dun "f^+f^ f 1Q f Q1 



7— ch 



'v— ch 



fmif: 
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f, 



(m+l)0 



fni 



(8) 



Using eq. 21 eq. [5] and eq. [7] we can express few factorial 
moments in terms efficiency and average of neutral pion 
fraction. 



fio 

foi 

f'll 

120 
102 



(1 - /) e ch (N) 
(f) 2e 7 (N) 

(/(l-/)) 2£ 7 e ch (N (N-l)) 
((l-f) 2 )e 2 h (N(N-l)) 
(f)ie 2 (N (N - 1)) + 2e 2 (/) (N) 



Substituting these in eq[T]we obtain 



7— ch 
dyn 



((i-/) 2 ) 



(f 2 ) 
177 



9 </(!-/)) 
L </)<!-/> 



(iV(JV-l)) 

W 2 



5I7P' 



(9) 



We note here that the for generic case the term inside 
the bracket is zero and wc have 



7— ch 
dyn 



l 



generic 2 (fj (-/V) 



(10) 



Using proper combination of factorial moments and do- 
ing a simple method of event mixing one can extract 
the generic value of v 1 ^^ (see appcndix- IX Al for details). 



Subtracting the generic value of v\ y ^ one can get rid of 
the last term in eq|U 



by 



So we propose a modified variable Vdyr. 



dyn 



given 



7— ch 
dyn 



'<(W) 2 > , if 2 ) .(/(i-/)) 



(N(N-l)) 
(N) 2 



(fY 



(/>(!-/) 



(11) 



In ideal scenarios when all the particles are detected 
one can approximate go(z.y) = zt and <7ch(z c h) = z c h- In 
that case one can show using cq|5] and cq[7] that 



7— ch 
dyn 



generic 2 (N) (f) y/ (N ch ) (Nj 



(12) 



irrespective of any value of Vdyn- So in that case the 
observable Avd yn can be estimated to be 



Av 



,7-ch 
dyn 



7— ch 
dyn 



(13) 



Following the same approach one can express the vari- 
able r TO ,i as 



7-ch = (/(1-/D(1-/) 
((1 -/)-+!> (/) ' 



(14) 



Now we would like to discuss the sensitivity of these two 
variables for a given fraction of DCC like signal. If re- 
fraction of events has DCC like domain formation, in 
simplistic case one can assume that the distribution of 
neutral pion fraction to be a combination of generic and 
DCC probability distribution given by 



1 



+ (l-x)5[f 



2\/7 

So for Ai/dyn we have from eqfTT 

((l-/) 2 ) (/ 2 ) 
(1 - ff </) 2 
(N(N-l)) 



(15) 



Ay 



,7-ch 
dyn 



(/(I-/)) 

(/>(!-/) 



{NY 
x (N(N 
V9 ~ 



1)> 



(Ny 



signal 



(16) 



which is proportional to the fraction of DCC-like events. 
Avdyn shows very high sensitivity to DCC like signal 
but it has dependency on the parent multiplicity and 
consequently to the collisions ccntrality. In later sec- 
tion wc would discuss this issue in detail. In case 
parent distribution is Poissonian, the fluctuation term 
(N(N - 1)) / (N) 2 would be equal to 1 giving AfJ~„ cl1 ~ 
i/(5/9). ' 

The robust observable expressed in eq03] would have 
a very particular x dependence given by 



7 — ch 



1 



—F(m,x) 



(m + 1) 

where the function F{m, x) is given by 
F(m, x) 



'2V m + 1 r(m+5/2) 
,3/ r(m+2) 



(17) 



(18) 



For ideal DCC case (x=l), the function F(m,x)=l for 
all values of m. That gives r m ,i — 1/(?7z-t~1). For generic 
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case(x = 0), r TO ,i=l. Fig.l shows the sensitivity of r m ^ 
for small signals of DCC. The functional form given in 
cq ll7l can be used to extract x from a fit of r m ,i with 
in. In the derivation of eg 1161 and cq |17l we have assumed 
that the parent multiplicity distribution are similar for 
both the generic and DCC case and the efficiency factors 
arc constant and independent of multiplicity and other 
kinematic parameters. 

III. EFFECT OF MIS-IDENTIFICATION 

There are additional complications in realistic scenar- 
ios that have not been taken care of in the above pre- 
scriptions. The study of 7 — ch correlation is often com- 
plicated by mis-identification of charge particles by pho- 
tons and vice versa. High energy depositions of charged 
hadrons can form a cluster in photon detector. Similarly 
photon conversion can show up as single or doubly de- 
tected tracks or clusters in charge particle detectors. In 
both the cases the observables are affected. Following the 
approach of the application of cluster decay theorem dis- 
cussed in previous section, we obtain the modified forms 
of the generating functions 

<7ch(£ch, Z 7 ) = (1 — £ch — £ch, 7 ) + Ech^ch + £ch,7^7 
5o(Zch, ;z 7 ) = ((1 — £ 7 — e 7jC h — £ 7 ,2ch) + £7X7 

~t~ ^7,ch Zch £7,2ch ^ch , (19) 



where we view neutral pions decay with 100% "efficiency" 
into two photons which themselves "decay" with a few 
modes. e c h and £ 7 are the efficiencies of detecting a 
charged particle and a photon, respectively. £ c h, 7 is the 
probability of charged particle being identified as a pho- 
ton cluster and £7,0111 £7,2011 are the probability of a pho- 
ton being identified as one and two charged particles, 
respectively. Substituting these in eqJS]one can calculate 
factorial moments folded with the contamination effect. 
The factorial moments are expressed as, 

fio = ((1 - /Kh + 2/ (e 7 , ch + 2£ 7 , 2 ch)) (AO 

foi = ((1 - /)£ch,7 + 2/£ 7 > (AO 

f u = (N(N - 1) ((1 - f)e ch + 2/ (e 7 , ch + 2£ 7 , 2ch )) 

X ((1 - /)£ ch , 7 + 2/£ 7 ) + 27V/£ 7 (£ 7 , ch + 2£ 7 ,2ch)> 

f 2 o = (N(N - 1) ((1 - f)e ch + 2/(£ 7)Ch + 2£ 7i2 ch)) 2 

+ 2Nf (2£ 7 , 2ch + (£ 7)Ch + 2£ 7 , 2ch ) 2 )) 
f 02 = (N(N - 1)((1 - /) £ch , 7 + 2/£ 7 ) 2 + 2iV/£ 7 ) (20) 

This would lead to very complicated (see appendix- IX Bp 
dependencies of Avd yn and r mj i on various efficiency fac- 
tors. However a relatively simple form can be obtained in 
the limit of small values of £ 7iC h and £ 7; 2ch- So in case of 
small photon conversion in the charged particle detector 
one can express Av c i yn as 



A 7— ch 
"dyn 



( 



\ 



((W) 2 ) 



((1-/^+2/ 



(1-/) ((l- 2/ 



( X - ^ ((1 - f)^ + 2/) 2 (1 - /) ((1 - fV-ff + 2/) 



(N(N-l)) 
(N) 2 



(21) 



where the the generic value of Vdyn in this case will be 
given by 



7— ch 



generic 2 (f) (N) ^ 6ch -~> -f lj 



(22) 



which is small number for large values of (N) . The robust 
variable r m l can be represented as 



7 — ch 



(a-/r ((i-/)^+2/)} (1-/) 

((!-/)-+!) ((i_ + 2/ 



(23) 



Unlike previous case it is not possible to eliminate the 
efficiency factors in eqEH and eql23l So in this case if 
we want to analyze the sensitivity of those variables to 
x-fraction of DCC signals, we can use eqlTSlto obtain the 



modified forms as 



Av 



7— ch 
dyn 



(N (N — 1)) 



r m ,i = 1 



5/9 j£o^ +l y (N) 
1 



mx 
m + 1 ( E ch,-y 



1 



F(m,x). (24) 



where F(m,x) is given by eqfTS] We can see that mis 
identification of charged particle as photon reduces 
the effective fraction of DCC events. The contami- 
nation factor appears as a ratio of £ c h.7/£7 keeping 
the functional form of the variables (cafTHl cqfTT]) un- 
changed. We note here that AvJ~^ h has quadratic 
dependence on contamination factor whereas ri.i is 
affected only by a linear factor. This is because Ai/J~^ h 
contains an extra photon fluctuation term absent in r m \. 
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(a)Variation of r m i with fraction of DCC signal 



(b)Variation of r m i with higher order moments 



FIG. 1: Sensitivity of variable r mi i and its higher moments. The higher order shows more sensitivity to small signals of 
anti-correlation. 



IV. RESONANCE EFFECT 

Resonance decays like p — > 7r ± 7 is equivalent to ar- 
tificial increase of pions and photons from generic case. 
Decays like lo —> ir° + tt^ 1 would give rise to correlation 
in the pions. Overall effect of resonance would be equiv- 
alent to event-by-event fluctuation of charged or neutral 
particles. The effect of resonance leading to increase in 
photon and charged particle multiplicity could be consid- 
ered to be equivalent to increase in efficiency of photon 
and charged particle detection. For event-by event fluc- 
tuations of efficiency would affect the observables, for e.g. 
the variable r TOj x given in eg 1141 will be modified as 



7 -ch _ (/(l-/) m )(W) (gTgS) fcc*) 



(1 - /) 



m+1 



<Cn +1 > <*r> 



(25) 



It is difficult to conclude the behavior of the variables 
from the above expressions without putting a realistic 
number for the efficiencies. To study the effect of reso- 
nances in a more detailed way (scc lVIIl IVIIII) we have 
used Monte-Carlo models in which resonances are in- 
cluded. 



CENTRALITY DEPENDENCE 



all the sources. If we assume heavy-ion collision to be 
a linear superposition of many identical nucleon-nucleon 
collisions, under identical source approximation we can 
calculate the centrality dependence of the variable using 
"central limit theorem" (CLT) [H. From CLT it follows 
that mean and variance of multiplicity would be given by 

/ Ns \ N s 

M(N) = M I ]T nA = ]T M(Ni) = N s M(Ni) 



N S 



a 2 (N)=a 2 \J2 N A = ^o 3 {N i )=N s a 2 (N i ).(26) 



So from CLT we have the dependence M(N) = aNs and 
<j(N) = P^/Ns- Let us assume N to be equal to the total 
number of produced pions where we have = aN c h + 
bNj. N could also refer to individual number of photons 
or charged particles. In that case similar argument also 
holds for M(N ch ) - ai N s and M{N 1 ) - a 2 N s . The 
variance of total numbers of pions would give 



(Nl) = ((aN ch + bN 7 ) 2 ) ~ fcNs + foN 2 s 



(27) 
(28) 



In heavy ion collisions, signals are expected to have 
centrality dependence, it is therefore important to study 
the centrality dependence of the 7 — ch correlation. In a 
heavy ion collision, let us consider Ns numbers of iden- 
tical sources are responsible for particle production. If 
Ni is the number of particles produced from i-th source, 
any variable V(iVj) will have a distribution identical for 



and if we express pion multiplicity in terms of charged 
and photons we get, 

(N 2 ) ^P^s + hN 2 
(N ch N 7 ) ex PtN s + p 2 N 2 s (29) 
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So from eqQ] and eq[2] we can calculate the centrality 
dependence of the observables. For v dyn one has 



7— ch 
dyn 



A 



B 



= A' 



B' 



(30) 



ch/ 



which is in fact the centrality dependence of all three 
terms in eqfTJ Here we note that the constants A' and B' 
could be cither positive or negative depending on which 
term in eqQ]is dominant. The variable Av dyn would have 
the similar centrality dependence which is obvious from 
the form of ea ll3l In heavy ion collisions, number of 
source participating in particle production can also be 
assumed to be proportional to number of participants 



part , 



of the collision. 



In that case ^X,„ ch is ex- 



pected to show a scaling behavior of the form A + B /X 
with X being either observed multiplicity or a Galuber 
variable N part . In case of experimental measurements 
it is more convenient to express fluctuation variables in 
terms of measured multiplicities. 

Based on similar approach one can comment on the 
centrality dependence of the robust observable. In the 
most general case one can have 



r m ,i 



_P 

m 



(31) 



which shows identical dependence in both numerator and 
denominator. So according to CLT, behavior of r mj i with 
multiplicity depends on the coefficients ai and Pi. How- 
ever it must be noted that breakdown of scaling from 
CLT would have several implications. The picture of 
identical source emission may not be valid in the case 
for formation of domains of DCC. In that case one might 
observe deviation from proposed scaling. 



VI. EFFECT OF MIXTURE OF PION SOURCES 

In this section wc would like to discuss the effect on 
the observables when event wise pion sources are inde- 
pendent of each other. So far we have considered that in 
a DCC event, all the pions detected in a given coverage 
are coming from the decay of the domains of DCC. This 
assumption might be valid when the detector coverage 
is same as the size of DCC domains. Let us consider 
a case when x-fraction of events analyzed has DCC like 
fluctuation carried by y-fraction of total pions. This is 
a realistic scenario when the size of the domain of DCC 
is smaller than the detector coverage. So for DCC pions 
we have (N) D = y (N) and for generic pions we have 
(N) G = (1 — y) (N), N being the total number of pions. 
The probability to find Njj pions carrying DCC signal 
will be given by P(N D , N, y) = n C Nd y N °{l- y) N 

d = y 2 (N(N-i)). 



which would give (N(N — 1)) 



Now 



in this case the generating function of cq|4] will be re- 
placed by 

Gobs = X' 'GdCC + %GvcgG g ener i c + (l—X — x')G 

generic 

(32) 

in which we view cases with 100% DCC production (x' 
fraction of events), 100% generic production and a mix- 
ture of two as three "decay modes" of a super clus- 
ter. Here Gdcc has the probability distribution V(f) = 
1/2a/7 and G generic has V(f) = S(f - 1/3). Since we 
think the case of 100% DCC production is the least real- 
istic, in the following we simplify our expression by tak- 
ing x' = 0. Now different factorial moments will become 
functions of x and y (see appcndix- IX CI for detail). In 
this case the observables are modified accordingly , for 
Av dyn from eq[T3] will be given by 



x 2 (N(N-1)) 



(33) 



which consistent with the expression eq |16l for y = 1 case. 
For Poisson like parent distribution A^yn can be ex- 
pressed as 



Av dyn = —y 2 . 



(34) 



We note here that Av dyn still shows the proportionality 
with the fraction of DCC events x. And the interesting 
fact is that quadratic dcpcndancc on y means Av dyn is 
more sensitive to the change of fraction of pions carrying 
DCC-like signals. 

In similar approach we can express n,i to be 



2xy 2 



xy' 



(35) 



This expression is consistent with the approximate ex- 
pression of r^i given in Ref. 24 1 for small values of x. The 
higher order moments will have corrections from higher 
orders of y which are smaller. To the lowest order ap- 
proximation, the expression given by ea !17l is still valid 
with fraction x replaced by xy 2 . 



7— ch 



1 



mxy 
(m + 1) 



F(m,xy ) 



(36) 



A functional fit of r m> \ with m to experimental data by 
the above expression can restrict the contours of x and 

y- 

VII. MODEL PREDICTION 

In this section we would like to study the behavior of 
observables from different models available to describe 
heavy ion data. There are theoretical predictions of 
isospin fluctuation for a statistical system of pions (25l . [26| . 
It can be shown that a system of Boltzmann gas of pions 
in the grand canonical ensemble (GCE), gives (N n o) = 
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FIG. 2: Multiplicity dependence of observables n t i and Aud yn as predicted from different models. The curves represent the 
results from different ensembles of Boltzmann gas of pions from eq |37l and eq |39l as described in the text. The markers are from 
different Monte-Carlo models. The error-bars are statistical. 



(N 7r ±) = z and one finds mean-square of pion multiplic- 
ity and charge-to-neutral pion correlation to be related 
to mean multiplicities as 

(iV^o) = (A^o) + (A^o> 2 

(N*±) = (N„±) + (N 7r ±) 2 

(iV ff oiV ff± ) = (N^)(N„±) (37) 

where z is the single particle partition function. In ref(26j 
it was shown that for an ideal scenario where one as- 
sumes the total isospin of the system to be zero, above 
mentioned relationships will become complicated. An en- 
semble of the total isospin 1=0 gives 

(N n o) = (N n± ) = j + j (38) 
but the mean-square pions multiplicities are modified as 



z 

15 



z 

10' 



(39) 



We can generalize this result and apply in case of our 
observables of 7 — ch correlation. The dependance on z 
can be eliminated and final observables can be expressed 
in terms of experimentally observed quantities like mea- 
sured multiplicity (say yj (iV c h) (N^)). In this case one 
has (iV 7 ) = 2(iV w o) and (JV ch ) = {N n+ + N % -) = 
2(7V 7I .±) . Also for decay of neutral pions, in case when 
all the photons are detected we have 2 er 2 w 2<7 2 „. One 



can express the mean-square multiplicity to be 

<7V 2 ) =4(7V 2 „) , (iV 2 h ) = 2(iV 2 ± ) + 2^+^-) 

(40) 

and the correlation term will be given by (iV 7 iV c h) = 
4 (N n oN n ±). Now we have 



120 
ho 

f()2 
f 2 " 

fn 



1 / (N„±(N„± -1)) | (jVjV) 



1 / (N^jN^ - I)) 



1 



(41) 



and eqUH we can estimate v 1 ch 



dy n 

Using eq[13] we 



So using eql3"Tl eq[ 
and ru for GCE and 1=0 systems, 
can estimate Av 1 ^ . For GCE we get from cq|37l 

and eqUU Vd yn = 1/ \J (N c h) (Nyj, which gives correct 
multiplicity dependence as predicted from CLT. So we 
for GCE. The system of 1=0 gives 



have Av 



,7-ch 
dyn 



AvJ y ™ 0.98/yJ(N ch ) (iV 7 ) which also agrees with 

the CLT predictions as shown in fig 2(b) In case of GCE 
n,i is predicted to be 2/(1 + 1 / yj (N ch ) (Ny)) which be- 
comes 1 for large values of multiplicity. For system of 
1=0, ~ 1 for all values of */ (iVch) (Ny) a s shown in 



fig 2(a 



For Poissonian case cr 7 = yj (Ny) = yjl (N^o) = v2cr^o 



We have estimated various observables and their cen- 
trality dependence using different monte-carlo eve nt g en- 
erators like HIJING [13, AMPT H| and UrQMD || for 
top RHIC energy. For our calculation we choose one unit 
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(b)Centrality dependance of variablc(markers) va yn 
and Avdyn and fits(lines) predicted from CLT. Here 
N c h and Nj refers to the mean multiplicities of 
charged particles and photons. 

FIG. 3: Prediction of variables from different models. 



of rapidity in forward direction 3 . We do the centrality 
selection based on putting cuts on impact parameter fol- 
lowing Glauber model calculation. FigJ5] shows the cen- 
trality dependence of the obscrvablcs. The variable r^i 
shows flat centrality dependence within error bars. As 
shown in fig 2(a) and fig |2(b)| the results from different 
monte-carlo models are consistent with each other and 
the values from the statistical model of Boltzman gas are 
consistent with other models towards higher multiplicity. 



3 both STAR and ALICE experiments has the setup of simultane- 
ous measurements of charged and photon in one unit of rapidity. 



At lower multiplicities they have qualitatively different 
nature probably due to presence of various other effects 
in the monte-carlo models. 

Fig 3(a) shows the variation of r TOj i with its order m. 



Results from all the model are consistent with the generic 
case of pion production. Fig |3(b)| shows the centrality 
dependance of Vdyn and Avd yn predicted from HIJING. 
For comparison of centrality dependance predicted from 
CLT, wc have fitted the points with functional form of 
A + B/y/(N ch ) (AT 7 ). This yields a value of A m 5 x 
10~ 5 and B = —0.6 for Ay^. We also note here that 
the sign of Ah>d yn is negative for low multiplicity. This 
means that the Raw HIJING includes some intrinsic 7 — 
ch correlation making the last term of eqQT]to dominate 
over individual fluctuation. This can be attributed to the 
resonance decays present in HIJING model. For DCC 
like signal sign of Avdyn should become positive for all 
centralities. 



VIII. DCC MODEL 

We have tried to implement DCC like anti-correlation 
signals in HIJING events. For a given event we changed 
the neutral pion fraction to follow I like distribution 
by flipping ir° to . And finally wc decay the neutral 
pions to photons. In the process of flipping wc make sure 
that the charge and isospin conservations are maintained. 
FigH] shows the /-distribution after the implementation 
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FIG. 4: Histograms showing distribution of neutral pion frac- 
tion for generic and DCC events from HIJING 

of DCC in HIJING. For generic event the neutral pion 
fraction is peaked at 1/3 and for DCC events it has a 
long tail. Since the variation of DCC like domain forma- 
tion with rapidity and azimuthal angle is not known, we 
perform this flipping for all the particles. This produces 
1/2 v 7 / like distribution over all phase space. To make 
the scenario more realistic we do the calculation of the 
final variables using total number of detected photons 
and charged particles rather than considering only pions. 
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1.05 




soft pions. The difference between the two models of 
string fragmentation and hydro is recently contested in 
ref.(3Q|. We therefore randomly choose pions produced 
in HIJING events, treat them to be thermal and im- 
plement l/2y/J distribution. Fig JS] shows the central- 



0.95 



0.85 



100 200 300 

V N c h xN Y 

(a) Variation of r±i with multiplicity 
0.4r 



400 



0.3 



0.2 



c 

>. 

73 0.1 

$ . 



-0.1 



-0.2 



o Raw HIJING 


AuAu 200 GeV 
2.5<r|<3.5 


— 














HIJING +DCC(%) 




0.1 

1 

5 






10 


100 200 


300 400 


V N ch 


X N y 



(b) Variation of i>dyn with multiplicity 

FIG. 5: Multiplicity dependence of observables n,i and Vdyn 
as predicted from DCC implemented HIJING model. Here 
JV c h and iV 7 denotes the mean multiplicities of charged parti- 
cles and photons for various centralities. The gray band shows 
the statistical error in model calculation. 
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(a)Value of r m i with m for various fraction of DCC 
events. The solid markers are when N c ^ and N y 
includes all the charged particles and photons and the 
hollow markers are when only pions are source of 

charged particles and photons. The curves are 
estimations from eg 1 171 and points are from DCC 
implemented HIJING. 
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Other dominant sources of photons and charged parti- 
cles include 77, charged kaons and protons respectively. 
It is difficult to extract the fraction of primordial pions 
on which the DCC-like probability distribution could be 
implemented. HIJING has minijet like environment in 
which the production mechanism are "string fragmenta- 
tion" and the abundance of particles are weighted by the 
spin giving large fraction of pions coming from decay of 
resonances. The primordial pions coming directly from 
string fragmentation are much smaller. Alternative envi- 
ronment like hydro models where the massive resonances 
are exponentially suppressed would give large fraction of 



I III I LJ l_l I I l_l LJ I i l_L 
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(b) Variation of r m . i in DCC events with fraction of 
pions coming from decay of DCC domains. Curves are 
estimations from cq |36l 

FIG. 6: Sensitivity of r m ,i to DCC like signals, estimation 
shown for x = 1. 

ity dependance of the two observables and their sensi- 
tivity for different fraction of DCC events, r^i shows 
almost flat dependance on multiplicity and we also find 
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similar depcndancc for all higher moments of r m ,i- Ab- 
solute values of rn are consistent with the prediction 
(n,i = (5— 2x)/(5+x)) from cqHTl For higher fraction of 
DCC the centrality dependance has slight non-monotonic 
behavior. This is also seen in Aud yn - As expected from 
cq llGl the values of Avdyn show proportionality with the 
fraction of DCC events. The absolute values of Avd yn 
are also very close to « x/(5/9) as predicted in calTrJl 
The centrality dependence causes ~ 15% variation of the 
values of most central to peripheral events for Avd yn - 
FigJBfa) shows that the variation of r m< \ with to. The 
results from model match the theoretical curve (cqfTT]) 
when one considers only pions as source of charged par- 
ticles and photons, however when all other sources are 
considered the results arc off towards lower side. A more 
detailed study of the sensitivity to fraction of DCC pions 
is shown in Figj6jb). We have shown the sensitivity of 
f m ,i with the fraction of detected pions carrying DCC- 
signals. In figEth) we also plot the curves obtained from 
cq 137)1 The effect of resonances present in HIJING seems 
to be resulting in reduced sensitivity of r TOj i for lower 
fraction of DCC pions. 



is proportional to xy 2 . For generic case of particle pro- 
duction from CLT, it is predicted to be inversely propor- 
tional to multiplicity. The sign of Avd yn would indicate 
the dominance of correlation over anti-correlation. 

We also discuss the applicability of the Minimax vari- 
able r mi i for heavy-ion collisions. r m i seem to be flat 
with centrality. Higher orders of r m ,i shows larger sen- 
sitivity x and can have contribution up to y m+l . A sim- 
plified form of the functional dependance of r OT) i with m 
has been calculated in generating function approach for 
lowest order of y 2 . This would be useful to restrict the 
signal strength xy 2 by fitting the experimental data. 
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IX. SUMMARY 



X. APPENDIX 



We have developed a procedure for generalization of 
methods for studying 7-charge correlation in heavy-ion 
collisions. One of the primary motivations for this study 
could be the search for DCC-like anti-correlation sig- 
nals relevant to the ongoing heavy ion program at RHIC 
and LHC. We have discussed the robustness of two vari- 
ables Avdyn and r m> i and have studied their central- 
ity(multiplicity) dependance. The sensitivity of the vari- 
ables have been studied with the fraction of DCC type 
events(a;) and the event wise fraction of DCC pions(y). 
These variables where evaluated from different models 
which do not include the physics of DCC. We have also 
developed a Monte-Carlo model where DCC domains 
have been implemented in HIJING to see the sensitivity 
of those variables with DCC signals. Our results show 
that the model predictions of the variables are consistent 
with the theoretical predictions using generating function 
approach. 

We have implemented the detector effects like effi- 
ciency of detection, mis-identificationto study the effect 
on observables. The mis-identification factor reduces the 
effective signal strength for which an approximate expres- 
sion has been derived in generating function approach. 
r m i has been found to be more robust towards mis- 
identification of photons as compared to Avd yn - The res- 
onance decay can induce correlation which can suppress 
the anti-correlating DCC signal. A quantitative idea of 
resonance can be obtained from Monte Carlo model that 
implements DCC in which we can vary the number of 
DCC candidates and see the sensitivity of r m ,i- 

We have seen that the variable Avd yn is highly sensi- 
tive to the fractions x and y. In a given centrality Avd yn 



Mixed events 



While analyzing data sample to calculate v 



,7-ch 
dyn 



can estimate the generic term by doing a mixed event 
analysis. A simple method we prescribe is to take to- 
tal number of photons and total number of charge par- 
ticles from different events would only effect correlation 
terms like fu in vj y ^ h keeping other factorial moments 
unchanged. In such case we must have 



fu 



fio foi 



(N(N-l)) 
(AO 2 



(42) 



Taking a particular combination of factorial moments we 
can calculate the generic value for Vd yn we need to cal- 
culate Avdyn For example one can show that 



fu 



. fa 



fio foi 



102 



P f 2 

I 10 M)l 



mixed 



2(f) (N) 



(43) 



which is equal to v 9 dy ^" c . But m case °f contamination 
effects present in the data sample one cannot apply this 
simple method since in that case the efficiency terms can- 
not be eliminated from Vd yn - A full GEANT simulation 
with a known event generator which doesn't include the 
physics of DCC is suggested to estimate the generic value 

Of Vdyn- 



B. Mis-identification 

In case of mis-identification of photon as charge parti- 
cles and vice-versa the fractorial moments are modified 
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as given in eqj20] The observables Avd yn and r m ,i will be given by 
I 



a 7 — ch 

dyn 



' fa I /)£ch + 2/£ ^ ch ) 2 ) | fa - fi 6 ^ + 2f£ ^) 2 (((1 - /)£ch + 2/ £7 , ch ) ((1 - /) £ch , 7 + 2/ £7 ))\ (N(N-l)) 
((1 - /Kh + 2/e 7 , ch ) 2 ((1 - /)e ch , 7 + 2/e 7 ) 2 ((1 - /)e c h + 2/e 7 , ch ) <(1 - /)e C h, 7 + 2/e 7 ) j (TV) 2 

(44) 

= (N(N - 1) ((1 - f)e ch + 2/e 7 , ch ) m ((1 - /) gc h, 7 + 2/e 7 ) + 2Nfs 1 e^ ch ) ((1 - f)e ch + 2/e 7 , ch ) 
(TV(TV - 1) ((1 - /) £ch + 2/£ 7 , ch ) m+1 + 2Nf (2e 7i2ch + e 2 ch ) ) ((1 - /)e ch , 7 + 2/e 7 ) 

(45) 



In case of £ 7 , c h = £ 7 ,2ch = one recovers eql2T1 and eql23l moments are given by 



C. Pion mixture 

In case of x-fraction of DCC events containing y- 
fractions of pions carrying DCC signal, different factorial 



f 10 = (l-/)e ch (iV) 
foi = (f)2e 7 {N) 

which is same as the case corresponding to y 
higher order moments are modified to be 



1. But 



fn = 
fao = 

fo2 = 



(2xy(l - y) (/) (1 - /) (TV) 2 + ((1 - *I/(2 - y)) (/(I - f)) G + xy 2 (f(l - f)) D ) (TV (TV - 1))) 2e 7 s 
(2xy(l - *,) (1 - ff (N) 2 + ((1 - xy(2 yj) ((1 - /)% + xy 2 ((1 - /)%) (TV (TV - 1))) £ 2 h 
y) (/) 2 (TV) 2 + ((1 - xy{2 - y)) (/% + xy 2 (f) D ) (TV (TV - 1))) 4 £ 2 + 2s 2 (/) (TV) 



2xy(l 



which gives 



2xy 2 



5 + xy A 

35 - xy 2 {2l-Ay) 



35 + xy 2 (2\ - 2y) 
and so on. The general formula for r m i is given by 

,2 



(47) 



r m ,i 



1 



mxy 
(m + 1) 



F(m,xy 2 ) + 0{xy A ) 



(46) 



r 



in which r TO i will have contribution up to xy m+ . Since 
y < 1 higher order contribution of y are smaller and the 
approximate form of the above expression would be given 
by 



mxy 
(m+ IT 



F(m, xy 2 ) 



(49) 



(48) where F(m, zy 2 ) is given by cqfTSl 
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